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Abstract An experiment that involves two distant mesoscopic SQUID rings is studied. 

The superconducting rings are irradiated with correlated photons, which are 
produced by a single microwave source. Classically correlated (separable) and 
quantum mechanically correlated (entangled) microwaves are considered, and 
their effect on the Josephson currents is quantified. It is shown that the currents 
tunnelling through the Josephson junctions in the distant rings, are correlated. 
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1. Introduction 

A fundamental property of superconducting quantum interference devices 
(SQUIDs) is that they exhibit quantum coherence at the macroscopic level [1]. 
This property may be used for the purposes of quantum information processing 
[2, 3]. 

A lot of research on superconducting devices investigates their interaction 
with classical microwaves. On the other hand the use of nonclassical microwaves 
makes the system fully quantum mechanical and interesting quantum phenom- 
ena arise. For example, in this paper we show that entangled two-mode mi- 
crowaves produce correlated currents in distant SQUID rings. 

Nonclassical electromagnetic fields at low temperatures (fceT ^ hu)) have 
been studied for more than twenty years both theoretically and experimentally 
[4]. The interaction of SQUID rings with nonclassical microwaves has been 
studied in the literature [5, 6]. 

In previous publications [7] we have studied the effects of entangled elec- 
tromagnetic fields on distant electron interference experiments. In this paper 
we review and extend further this work in the context of SQUID rings. We 



consider two mesoscopic SQUID rings, whicli are far from each other and are 
irradiated with entangled microwaves, produced by a single source (Fig. 1). 
It is shown that the Josephson currents in the distant SQUID rings are corre- 
lated. The photon-induced correlations between the currents are quantified. It 
is shown that the current correlations depend on whether the photons are clas- 
sically correlated (separable) or quantum mechanically correlated (entangled). 
The difference between separable and entangled microwave density matrices [8] 
is in the nondiagonal elements; and the effect of these nondiagonal elements on 
the Josephson currents is explicitly calculated. 

2. Interaction of a single SQUID ring with nonclassical 
microwaves 

In this section we consider a single SQUID ring and study its interaction with 
both classical and nonclassical microwaves. 

For irradiation with classical microwaves, the Josephson current is I a ~ 
7i sin^A, where ^a = 2e$A is the phase difference across the junction due to 
the total flux <i>A through the ring. We assume the external field approximation, 
where the back reaction (the additional flux induced by the SQUID ring current) 
is neglected; i.e., the flux CI a, where C is the self-inductance of the ring, is 
negligible in comparison to <1?a. The magnetic flux has a linear and a sinusoidal 
component: 

$A = VAt + <^A; (j)A= AsUl{LUxt). (1) 

Consequently the observed current is 

Ja = /i sin[ciJAi + 2eAsin(a;if)]; u)A = '2eVA- (2) 

We now consider the interaction of a SQUID ring with nonclassical microwaves, 
that are carefully prepared in a particular quantum state and are described by a 
density matrix p. The dual quantum variables of the nonclassical field are the 
vector potential Ai and the electric field Ei. Integrating these over the SQUID 
ring we obtain the magnetic flux and the electromotive force operators = 
§^ Aidxi, Vemf = §^. Eidxi. 

In the external field approximation the flux operator evolves as 

0(f) = ^2^ ' [a cxp(jcijf) + fflcxp(— ia;t)], (3) 

where ^ is a parameter proportional to the area of the SQUID ring and the a^, a 
are the photon creation and annihilation operators. Consequently the phase dif- 
ference 6a is the operator 

^a = <j-'Af + g[a exp(ia;i) + aexp(— jcjf)], g = v2e^; (4) 

and the current also becomes an operator, I a ~ Ii sin 9 a- Expectation values 
of the current are calculated by taking its trace with respect to the density matrix 
p, which describes the nonclassical electromagnetic fields, 

(/a) = Tr(pJA) = /iIm[exp(ia;Af)Vy(AA)], (5) 

Aa ~ iq exp{iujit) . (6) 
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W{x) is the Weyl function [9] given by 



W{x) = Ti\pD{x) 



D{x) — exp(a;a — x*a) 



(V) 



where D{x) is the displacement operator. Higher moments of the Josephson 
current quantify the quantum statistics of the electron pairs tunnelling through 
the junction. 





Figure 1. Two distant mesoscopic SQUID rings A and B are irradiated with nonclassical 
microwaves of frequencies uj\ and (j;2, correspondingly. The microwaves are produced by the 
source Sem and are correlated. Classical magnetic fluxes VXt and Vet are also threading the 
two rings A and B, correspondingly. 



3. Interaction of two distant SQUID rings with entangled 
microwaves 

In this section we consider two mesoscopic SQUID rings far apart from each 
other, which we refer to as A and B (Fig. 1). They are irradiated with correlated 
microwaves. Let p be the density matrix of the microwaves, and pA = Trep, 
Pb = TiAp, the density matrices of the microwaves interacting with the two 
SQUID rings A, B, correspondingly. When the density matrix p is factorizable 
as pfact = Pa ® Pb the two modes are not correlated. If it can be written as 
Pscp = X/ Pip^i®pBi, where pi are probabilities, it is called separable and the 
two modes are classically correlated. Density matrices which cannot be written 
in one of these two forms are entangled (quantum mechanically correlated) [8]. 

The currents in the two SQUIDs are 



(/a) = /iTr(pA sin 61a), (/b) = /2Tr(pB sin 61b). 



(8) 



The expectation value of the product of the two current operators is given by: 
(/a/b) = /i/2Tr(psin0AsineB). (9) 

In general (IaIb) is different from {I a) (Ib) and we calculate the ratio 

(/a/b) 



R^ 



{Ia){Ib} 



(10) 



For factorizable density matrices pfact = PA® PB we easily see that -Rfact = 1- 
For separable density matrices pscp the ratio -Rscp is not necessarily equal to one 
and numerical results for various examples are shown below. 
We also calculate the second moments 

(/i) = J?Tr[pA(sineA)2], (Jl) = j|Tr[pB(sinfe)']. (11) 

The statistics of the photons affects the statistics of the tunnelling electron pairs, 
which is quantified with the (Ja/b), {IJ^, (Ib) (and also with the higher mo- 
ments). 

3.1 Microwaves in number states 

We consider microwaves in the separable (mixed) state 

Pscp = i(liViiV2){iViiV2| + |iV2iVi)(iV2iVi|), (12) 

where A''i 7^ N2. We also consider microwaves in the entangled state \s) = 
2~^ "(|iViiV2) + \N2N1)), which is a pure state. The density matrix of \s) is 

Pent = P.cp + ^{\NlN2){N2Nl\ + \N2N1) {NlN2\) , (13) 

where the psop is given by Eq. (12). It is seen that the pont and the psop differ 
only by the above nondiagonal elements. 

In this example, the reduced density matrices are the same for both the sepa- 
rable and entangled states: 

Pscp.A = pent.A = Psop.B = Pont.B = ^(|A''l)(A^l| + \N2){N2\). (14) 

Consequently in this example (/A)sep = {/A)cnt, and also (Jb)scp = (/B)cnt. 
For the density matrix pscp of Eq. (12) we find 

(/a) = y exp f --^ j [LiVi(g^) + L]V2(q^)]sin(a;At), (15) 

(Ib) = y cxp f --2- j [iiVi(g^) + L]V2(g^)]sin(a;Bt), (16) 

where the L"{x) are Laguerre polynomials. The currents (/a), (^b) are in this 
example independent of the microwave frequencies uji,u>2. 

The expectation value of the product of the two currents [Eq. (9)] is 

{/a/b)scp = Iil2 0xp{—q )LNi{q )LN2{q )sm{ujAt)sm{LJBt). (17) 

Consequently the ratio R of Eq. (10) is 

p 4Ljv,(g^)LAf3(g^) ,,„. 

In this example the Rscp is time-independent. 
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The moments of the currents, defined by Eq. (11), are also calculated: 

{il) = ^{l"icxp(-2g2)[Ljv,(4g2)+Ljv,(4Q=^)]cos(2c^At)}, (19) 

{/!) = :|{l-icxp(-2g2)[iA.,(4g2) + LA.,(4g2)]cos(2u;Bi)}. (20) 

For the case of pcnt the (/a), (Jb) are the same as in Eqs. (15), (16); and the 
(Ja), {I'bI are the same as in Eqs. (19), (20). However the {IaIb) is 

{lAl'B)cnt = {IaI'b)scp + /cross, (21) 

where 

/cross = — exp(-g^)L^J"^i(g^)L^J"^^(g^)[cos(c<JAt + i.<jBt) 

_(_;^)iVi-iV2 gog(^^i _ ^gi)] cos(ni), (22) 

n = {Ni-N2){LOl-i02). (23) 

It is seen that the effect of entangled microwaves on Josephson currents is dif- 
ferent from the effect of separable microwaves. In this case the ratio R of Eq. 
(10) is 

-Rent = -Rsep + "°°^-; , (24) 

(/a)(/b) 
which is a time-dependent quantity oscillating around the _Rscp . 

3.2 Microwaves in coherent states 

We consider microwaves in the classically correlated state 

P-P = ^{\AiA2){AiA2\ + \A2Ai){A2Ai\), (25) 

where \Ai), \A2) are microwave coherent states. We also consider the entangled 
state \u) =J^{\AiA2) + | ^2^1)), with density matrix 

Pent = 2AAVscp +Af'\\AlA2){A2Al\ + \A2Aj){AlA2\), (26) 

where the normalization constant is given by 

AA= [2 + 2exp(-lAi-A2l')]"'^'. (27) 

For microwaves in the separable state of Eq. (25) the reduced density matri- 
ces are 

pscp.A = Psep,B = -{\Ai){Ai\ + \A2){A2\), (28) 

and hence the current in A is 

(/a)sop = -— exp(— -— ){sin[c<JAi + 2g|^i| cos(LJii — 6li)] 

+ sin[cJAt + 2gl^2lcos(uJii-6'2)]}. (29) 
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Figure 2. Rscp against (cji — iJ2)t for the number state of Eq. (12) with A^i = 1, Ai'2 = 4 
(line of circles), and the coherent state of Eq. (25) with Ai = 1,A2 = 2 (solid line). The 
photon frequencies are oji = 1.2 x 10^'' and L02 = 10^'*, in units where ks = h = c = 1. 



where Oi — arg(Ai), and 82 — arg(^2). A similar expression yields the 
current in B. We have also calculated numerically the ratio i?scp- 

For microwaves in the entangled state of Eq. (26) the reduced density matri- 
ces are 



Pont.A = Pcnt,B = J^\\Al){Al\ + \A2}{A2\ + T\A-i){A2\ + T* \A2} {Ai]) , 

where r = (Aij^a) = oxp(-|Ai|V2 - i^2JV2 + AIA2). 
The current in A is 



(30) 



(/A)cnt = 2N\lA)..p+M'EFi CXp ( -:|- ) h, 



(31) 



where £; = exp[-l^ip - \A2\'^ + 2\AiA2\cos{ei - 6'2)], and 

Fi = [exp((jj^i|S'A,i -q\A2\SA,2) +oxp(-g|Ai|S'A,i 
+q\A2\SA,2)] shiiuJAt + g| Ai|Ca,i + g| A2ICA.2), 



(32) 



with Sa,i = sin((j;it — 6*1), Sa,2 ~ sva{u]it — 62), Ca,i = cos(tJii — 61), 
Ca,2 = cosiujit — 62). A similar expression yields the current in B, and we 
have also calculated numerically the ratio Rent ■ 



3.3 Numerical results 

In the numerical results of Figs. 2 and 3 the microwave frequencies are ui = 
1.2 X 10^'^ ,u!2 — 10^'', in units where ks ~ h = c — 1. The critical currents 
are Ii — I2 ~ 1. The other parameters are ^ = 1, cja — Lui,iL>B = u:2, Ni = 1, 
N2 = 4, and 9i = 62 ~ 0. For a meaningful comparison between microwaves 
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(a)A=1. A =2 




Figure 3. (a) (Ja)scp - (/A)cnt, and (b) (iDsop - (^Dont against {iOi-L02)t for the coherent 
state Psep.A of Eq. (28) and pcnt.A of Eq. (30) with Ai — 1,A2 — 2. The photon frequencies 
are wi = 1.2 x 10^'' and L02 = 10^'', in units where ks = h = c = 1. 



in number states and microwaves in coherent states, we take them to have the 
same average number of photons, l^i]'^ — Ni and j A2P = N2. 

In Fig. 2 we plot -Rscp against {iLii—u>2)t for currents induced by microwaves 
in the number state of Eq. (12) with A^i = 1, N2 — 4 (hne of circles), and the 
coherent state of Eq. (25) with Ai = 1,A2 — 2 (solid line). It is seen that two 
different separable photon states with the same average number of photons give 
rise to different correlations between the SQUID currents. 

In Fig. 3 we plot (a) {/a)scp - {/A)ent, and (b) {I'Dscp - (/Dent, against 
(wi — uJ2)t for microwaves in the coherent state pA,sop of Eq. (28) and pA,ont 
of Eq. (30) with Ai = 1, A2 = 2. In Fig. 3(a) it is seen that the losephson 
current in SQUID ring A is different for irradiation with separable and entan- 
gled microwaves in coherent states. In Fig. 3(b) it is seen that irradiation with 
separable and entangled coherent states leads to different second moments of the 
current, which implies that the quantum statistics of electron pairs tunnelling the 
Josephson junction of SQUID ring A are different in these two cases. 



4. Discussion 

We have considered the interaction of two distant SQUID rings A and B with 
two-mode nonclassical microwaves, which are produced by the same source. 
The flux, the phase difference and the Josephson currents are operators and their 
expectation values with the density matrix of the nonclassical microwaves give 
the physically observed quantities. We have assumed the external field approxi- 
mation, where the electromagnetic field created by the Josephson currents (back 
reaction) is neglected and we have calculated various quantities. 

It has been shown that the Josephson currents in the two rings are corre- 
lated in the sense that (IaIb) is different from (7a) (/b) (for non-factorizable 



density matrices). We have considered examples where the photons are classi- 
cally correlated and quantum mechanically correlated; and we have shown that 
the non-diagonal terms in the latter case affect the Josephson currents. Further 
work in this direction could be the formulation of Bell-type inequalities for the 
Josephson currents, which are obeyed in the case of separable microwaves and 
violated in the case of entangled microwaves. 
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